A circumstellar disc in a high-mass star forming region 
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ABSTRACT 

We present an edge-on Keplerian disc model to explain the main component 
of the 12.2 and 6.7 GHz methanol maser emission detected toward NGC7538- 
IRS1N. The brightness distribution and spectrum of the line of bright masers 
are successfully modeled with high amplification of background radio continuum 
emission along velocity coherent paths through a maser disc. The bend seen in 
the position- velocity diagram is a characteristic signature of differentially rotating 
discs. For a central mass of 3OM , suggested by other observations, our model 
fixes the masing disc to have inner and outer radii of ~ 350 AU and ~ 1000 AU. 

Subject headings: star formation - masers - Interstellar medium - discs 

1. Introduction 

Disks are expected to form during protostellar collapse, and low-mass stars seem to 
provide good observational evidence for the existence of disks (e.g. Qi et al. 2003; Fuente 
et al. 2003). The situation is less clear for high- mass stars. While in several cases velocity 
gradients in massive star forming regions have been detected on large scales (> 10, 000 AU, 
e.g. Sandel et al. 2003), evidence for compact discs on scales < 1000 AU remains sparse. 
One possible example is IRAS20 126+4104, where there is good evidence for disc-outflow 
geometry around a 24M Q central object (Cesaroni et al. 1999; Molinari et al. 2000). The 
case for a disc around G 192.16-3.82, a protostar of ~15M (Shepherd et al. 2001), is less 
compelling. 

Class II methanol maser emission, a signpost of high mass star formation (Minier et al. 
2001), offers a potential indicator of discs since it often shows linear structures in both space 
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and position — line-of-sight velocity diagrams (Norris et al. 1998; Minier et al. 2000). These 
masers have been modeled as occurring at fixed radii within edge-on discs, although it is 
still unclear whether they arise primarily in discs or outflows (DeBuizer 2003; Moscadelli 
et al. 2002). One of the most striking examples of a maser line in both space and velocity is 
found in NGC7538-IRS1 N (Minier et al. 1998, 2000). We present here the first quantitative 
Keplerian disc analysis of this maser without invoking the assumption of a single radius. We 
find compelling evidence for the disc interpretation in this case. 
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Fig. 1. — Left: Zero moment (integrated emission over velocity) maps of the 6.7 and 12.2 GHz 
methanol masers in the main spectral feature of NGC7538-IRS1 N . Right: (8, v )-diagrams of 
the two masers, in which the data have been averaged over the spatial direction perpendicular 
to the major axis position angle. 



2. Maser data and other observations 

The NGC7538-IRS1 star formation region has been extensively studied at several wave- 
lengths. The spectral type of the central star has been estimated as 06 (Willner 1976; Camp- 
bell & Thompson 1984), in agreement with the SED from 2 mm (Akabane et al. 2001) to 
850 fim (Momose et al. 2001) which implies a central object with luminosity of > 8.3 x 10 4 L Q 
and mass of about 3OM . VLA observations uncovered an associated ultra-compact HII re- 
gion with peak brightness temperature 10,000-15,000 K (Campbell 1984; Gaume et al. 1995). 
Radio free-free continuum and recombination line observations show a N-S elongated struc- 
ture attributed to a large scale outflow. 
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The first VLBI observations of the 12.2 GHz methanol masers (Minier et al. 1998, 2000) 
revealed several groups of spots corresponding to different features in the single dish spectra. 
The subject of our study is the maser in the main spectral feature between about —56.50 
and —56.82 kms -1 . The 12.2 GHz data cube from published data by Minier et al. (2001), 
made from 1998 VLB A 1 observations, is shown in Figure 1. The interpretation of the linear 
structure in space and velocity as an edge-on disc (Minier et al. 1998) is supported by the 
positional agreement of the maser line with the waist of the hour-glass shaped continuum 
radio emission (Gaume et al. 1995) and a band of free-free absorption. The same structure 
is evident in the 6.7 GHz maser data, obtained in February 2001 with the European VLBI 
Network (EVN 2 ), shown in the figure. Taking into account the different angular resolution, 
the velocity-integrated emission of the two masers is similar out to ±20 mas in R.A. (Figure 
1, left). In addition, the 6.7 GHz data shows two disconnected 'outlier' maser features which 
lie close to the extrapolation of the maser line. The right panels in Figure 1 show the 
position-velocity diagrams for the two masers. They both show the same velocity gradient 
near the center and a bend at ~15 mas. 

The 12.2 GHz line involves the 2o — > 3-±E transition, the 6.7 GHz line the 5i — > 6qA + 
transition. Since E- and A-methanol can be regarded as distinct molecules for all practical 
purposes, the two masers can be expected to probe different physical conditions even when 
they reside in the same region. Yet in spite of these differences, the brightness peaks of the 
two lines are perfectly coincident in both velocity and spatial location. The agreement is 
within the 12.2 GHz angular resolution of 2-3 mas and the spectral resolutions of 0.048 and 
0.09 km s _1 p er channel for the 12.2- and 6.7 GHz data, respectively. Convolving the 12.2 GHz 
data to the 6.7 GHz spatial resolution we find that the ratio of the two emissions is constant 
to within 10% over the entire feature. 



3. Modeling 

The remarkable coincidence in space and frequency of the two masers is in fact a natural 
consequence of amplification of a background source by an edge-on rotating disk. In that 
case the maser brightness at displacements 9 and v from the position and velocity of the 
line-of-sight (LOS) diameter is 1(9, v) = Ib e T ^ ,v \ where Ib is the background continuum and 
r(9, v) the maser (negative) optical depth. Irrespective of the disk annuli where the masers 
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reside, the peak emission will always occur at (9,v) = (0,0), which picks out the largest 
optical depth To- With brightness temperature of 10 4 K for the background continuum, the 
function t(9,v) = In [1(9, v)/Ib] is a measured quantity. We find r = 18.32 and 15.99 for 
the 6.7- and 12.2 GHz masers, respectively. 

Consider an edge-on rotating disc at distance D and a point at radius p — r/D along 
a path with displacement 9. The rotation velocity V(p) and its LOS component v obey 
v 1 9 = V '/ p = f2, the angular velocity. In Keplerian rotation oc p~ 3//2 with f2 Q = fi(p ) = 
Dy^GM/R^ where p Q = R /D is the outer radius and M the central mass. Assume a 
Gaussian frequency profile for the maser absorption coefficient and denote the radial variation 
of its magnitude at line center by the normalized profile r](p) (J 7] dp = 1). Then 



as constant, although it could vary with p because of temperature variation and maser 
saturation. Saturation typically requires r > 10-15 across the disc radius, and the broadening 
of the absorption profile is then proportional to length in excess of this threshold (Elitzur 
1992). If the NGC7538-IRSlNdisc is located in front of the background source then the 
overall amplification with r ~ 18 implies an optical depth of only 9 across the radius, too 
low for saturation. If the radio continuum is centered on the star then the optical depth 
across the radius is ~ 18 and saturation can be expected to have an effect, but only at the 
very outer segments of the disc where a decrease in temperature could have an opposite, 
offsetting effect. We will consider these effects in a future publication. 

We parameterize 77 with a power law p~ p from an inner radius p x = R\/D to p D . Both tq 
and At>D are set directly from r(0, v). Modeling the full r(9, v) requires two free parameters 
for the integrand (Q and p) and two for the integration limits (p Q and h = pi/p Q ). Figure 2 
shows contour plots of t(9, v) for a representative Keplerian disc around a 30 M Q protostar 
with the At>D determined from the 12.2 GHz data at 9 = 0. The system is viewed edge-on 
at the adopted distance to NGC7538. 

The lowest contour (3%) outlines the angular velocity Q = v/9 on the disk boundaries. 
The shallow segments of this contour (lower branch for 9 < 0) trace out the disk outer radius. 
The steeper ones trace the inner radius when \9\ < p\ and the mid-line when \9\ > p u where 
the Keplerian v oc l^l -1 / 2 is evident. Along any path at \9\ < p h the velocity varies from 
t> min = Q Q \9\ on R to f max = Qi\9\ on Ri. As long as t> max — t> min < Af D , the full path 
remains velocity coherent and the maximal r is centered on the velocity determined from 
dr(v , 9)/dv = 0. This yields v = Cl9 where f2 = J ilr/dp, an average dominated by the high 
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Fig. 2. — Contour plots of t(9,v) (see eq. 1) for a Keplerian disc with Av-£,= 0.40 kms 1 
around a 30M Q star, viewed edge-on from D = 2.7 kpc. Constant maser absorption coefficient 
between radii R { = 378 AU (0714) and R Q = 1080 AU (074). 



angular velocity of the inner regions; for example, constant r\ gives = 2VL- X h/(1 + h 1 ^ 2 ). 
The linear relation between v and 9 is evident in the uniform inclination of the innermost 
contours. 

The condition t> max — t> min = At> D is met at a displacement 9^, and r gets contributions 
from only a fraction of the path when 9 > 9^. Since the longest coherence is at the outer 
radius, where Q is smallest, the r-contours bend toward lower velocities as is evident in 
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the figure around 9 k 0'/02. This bending reflects the change from inner-edge to outer- 
edge dominance of longest coherence. As 9 increases further beyond p\, the longest coherent 
segment moves from the outer radius to the disc mid-line. The maximal r begins to increase, 
creating the local outer maximum evident in the 50% contours at 9 ~ (X'39 before final 
truncation at 9 = p a . 

This discussion shows that the contours of highest amplification, down to 70% of To, 
contain dependence on p, h and Q but not on p Q . This has two important consequences: (1) 
Successful modeling of the strongest maser emission requires one less parameter. (2) Since Q 
involves only the combination M/R 3 , the central mass cannot be determined without obser- 
vations of (exponentially weaker) radiation close to the phase-space boundary to determine 

Po- 

Brightness measurements with dynamic range / map the r-distribution between r and 
r -In/. With r = 15.99 and / ~ 100, the 12.2 GHz data traces the T-contours in NGC7538- 
IRSlNdown to 70% of peak. Therefore the data shown in the lower left map of figure 3 
cover only the two innermost contours in figure 2, and if our model is applicable it should 
require only the three parameters h, Q and p. In fact, the first two can be determined 
quite accurately even without detailed modeling since the bend in the T-contours occurs at 
v k ~ ±Av D (1 + /i 3 / 2 )/(l - h 3 / 2 ) and 9 k ~ tt m v k , where fi m = + fi D ). We find directly 
that h = 0.39 and fi m = 0.037 km s^/mas, implying M '/ ' R\ = 3OM /(1OOOAU) 3 and leaving 
p as the only free parameter to match all other data points in the (9, v )-plane. Sensible fits 
to the 12.2 GHz data are obtained when 0.035 < Q m < 0.045 and 0.3 < h < 0.4. 

The lower right panel of figure 3 shows the model map for p = —0.5, the top panels 
show how well the fitting of t(9, v) reproduces the brightness profile and the spectrum. The 
negative p means that rj, i.e., the maser absorption coefficient, increases with radius. The 
relation between 77 and overall density involves the product of the methanol abundance, the 
fraction of methanol in the maser system and the fractional inversion. Any of these factors 
might give more weight to the disk outer regions. Considering the simple power law with 
sharp cutoffs we employ, a 70% increase in r] across the disk is compatible with a constant 
density profile. 

Our model is equally successful in fitting the main 6.7 GHz maser emission, as is ev- 
ident from figure 4. The value of p is the same for this maser, the only slightly different 
parameters are t = 18.32, /i=0.38 and tt m = 0.037 km/s/mas. The value for At> D was taken 
from the 12.2 GHz maser, which is considered to be more reliable because of the higher 
spatial resolution. Good fits of the 6.7 GHz data are obtained when 0.035 < fl m < 0.045 and 
0.35 <h< 0.45. The virtual identity of the disc segments of such different masing transitions 
is an important challenge for class II methanol pumping calculations. 



-7- 



4. Discussion 

The analysis presented here shows that a Keplerian disk explains successfully, down to 
minute details, the main methanol maser feature in NGC7538-IRS1 N . The angular velocity 
we find implies a disk outer radius of ~1000AU if the central mass is 30M Q , as determined 
by Campbell (1984). The methanol maser emission covers a substantial area of the disk, 
with the radius varying by a factor of 3. 

The t(9, v )-contours display a bend in the position- velocity diagram. Since v — 9, a 
straight line in this diagram implies a constant fl(p), i.e., either a single radius or, in the 
case of an extended region, solid-body rotation. Our analysis therefore uncovers a unique 
signature of differential rotation in the NGC7538-IRS1 N disc. The emission's smooth struc- 
ture and the observation's high sensitivity are crucial for this detection, which is made well 
inside the ~ 1.2 kms" 1 width of the spectral feature (see top right panels in figures 3 and 
4). These requirements preclude detection of the effect in extragalactic sources and explain 
why it could not be considered in the analysis of the H 2 maser disc in NGC4258 (e.g., 
Watson & Wallin 1994). However, a bend in the position- velocity diagram has been found 
in thermal line emission from low-mass stars and interpreted as steep-gradient solid-body 
rotation inside a shallow-gradient Keplerian outer region (Beckwith & Sargent 1993; Ohashi 
et al. 1997; Belloche et al. 2002). Since maser and optically thin thermal emission, which 
is directly proportional to optical depth, share the same kinematic structure, our analysis 
shows there is no need to invoke additional components on top of Keplerian rotation. 

Our analysis is valid for a radio continuum source placed either behind the maser disc 
or at its center. The latter possibility would fit well in the context of disc photoevaporation 
models (Hollenbach et al. 1994). In that case, the radial extent of disc photoevaporation by 
the central star in NGC7538 covers the entire maser region and could activate both methanol 
formation and the observed outflow, responsible for the radio continuum. The maser location 
at the outflow waistline makes this scenario especially attractive. 

The methanol maser data determines detailed disc properties with high precision yet 
the central mass remains undetermined since it requires the outer radius. And p Q is still 
unobservable because the angular extent of the detected emission does not trace the full 
disk, instead it is controlled by velocity coherence and the dynamic range of the observations. 
The amplification opacity at the disk outer edge (see fig. 2) is 50% of the central peak level. 
With r = 16 (18) at the center, detection of the outer peak requires a dynamic range of 
3000 (8100) if the radio source is located behind the maser disc, placing it beyond the reach 
of the current observations by a factor of ~ 30 (80). The requirements become even more 
severe if the continuum is centered on the star. The 6.7 GHz outlier features therefore cannot 
correspond to the outer peak of our model disc. Instead, these outliers are probably the first 
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detection of enhanced amplification by small inhomogeneities in an otherwise rather smooth 
disc. Such non-axisymmetric distortions of the disk surface, and their possible effect on the 
pumping, have been proposed by Durisen et al. (200f). It should be possible, though, to 
determine p Q by probing the boundary of the (6, v) region at 9 < p along the f-axis or 
another direction such as the normal to fl. We will study these issues in detail in a future 
publication. 

This study provides one of the few pieces of evidence to date for a compact (< 1000 
AU) disc around a massive protostar. It seems certain, though, that not all class II methanol 
masers occur in discs. Even in NGC7538-IRS1 N itself, the weaker methanol masers to the 
south, which are associated with other maser species (OH, H2O etc), are likely to arise from 
an outflow. The interpretation of ordered lines of methanol masers in terms of disc or outflow 
models must be decided case by case with detailed modeling. 

M. Pestalozzi thanks his fellow PhD student Rodrigo Parra for the useful discussions and 
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Fig. 3. — Modeling the 12.2 GHz data with a disc with the indicated parameters and Ai>r>= 
0.4 kms' 1 . Top: 6- variation of frequency-integrated brightness (left) and v- variation of the 
flux (right). In both diagrams the log of the quantities is shown. Solid lines show the data, 
dashed lines the model. Bottom: Contour plots of 1(6, v), the data to the left, the model 
to the right. The contours of the model have been obtained from the exponential of eq. 1. 
Contours are 5,10,30,50,70,90% of peak. The slope of the diagonal line is f^in kms _1 / mas - 
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Same as figure 3 for the 6.7 GHz maser. 



